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We derive a hierarchy of continuous-variable multipartite entanglement conditions in terms of 
second-order moments of position and momentum operators that generalizes existing criteria. Each 
condition corresponds to a convex optimization problem which, given the covariance matrix of the 
state, can be numerically solved in a straightforward way. The conditions are independent of partial 
transposition and thus are also able to detect bound entangled states. Our approach can be used 
to obtain an analytical condition for genuine multipartite entanglement. We demonstrate that even 
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optimization it is also possible to numerically verify genuine entanglement of some experimentally 
realizable states. 
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I. INTRODUCTION 

In the multipartite case there are many different no¬ 
tions of entanglement, ranging from the most specific 
to the most general ones. A specific kind of entangle¬ 
ment means that we precisely specify the groups of par¬ 
ties that are separable from each other, i.e., we specify 
a partition of the set of indices {1,... ,n}. A partition 
I = {/i,..., Ife} is a disjoint set of nonempty subsets of 
the indices whose union is equal to {1,..., n}. For exam¬ 
ple, in the case n = 4 there are partitions {2} U {1,3,4}, 
{1} U {2} U {3,4} and many others. We will use a more 
compact notation for partitions and write them as 2|134 
and 1|2|34. A partition I' is finer than a partition I (or 
the partition I is coarser than I') if for any J' there is an 
Ij such that /' C Ij. For example, the partition 1|2|34 
is finer than the partition 12|34. There are two extreme 
partitions, the trivial partition 12 ... n and the partition 
1|2|... |n. The former corresponds to the case where no 
information about separability properties is known and 
the latter corresponds to the notion of full separability, 
where all parties are separable from each other. The set 
of partitions with the finer-than relation is referred to 
as partition lattice and can be visualized as the graph 
shown in Fig. for n = 4. At the bottom is the trivial 
partition 1... n, then on the second line are all 2"“^ — 1 
partitions into two parts, next are all partitions into three 
parts, and so on. At the top is the partition 1|2| .. .\n. 
Partitions with k = 2 play a special role and are called 
bipartitions. 

From these specific kinds of separability one can con¬ 
struct more general types by considering mixtures of 
specific kinds according to some criteria. For example, 
a multipartite state is called fc-separable, k > 2, if it 
is a mixture of states each of which has k separable 
groups of modes (the corresponding partitions are on 
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FIG. 1. (Color online). The partition lattice of four parties. 
Nodes are partitions and edges correspond to the ’’finer than” 
relation. 

the same line in the lattice diagram). In the case of 
n = 4, a 3-separable state is a state which is a mixture 
of I|2|34-, 1|3|24-, I|4|23-, 2|3|14-, 13|2|4- and 12|3|4- 
separable states. The notion of n-separability of an n- 
partite state is the same as the notion of full separability. 
A 2-separable state is referred to as biseparable and the 
notion of biseparability is the most general of all — a 
state which is separable in any sense discussed above is 
automatically biseparable. On the other hand, full sepa¬ 
rability is the most specific kind of separability — a fully 
separable state is separable in any other sense. 

In general, a biseparable state is a mixture of 2”“^ — 1 
states corresponding to all nontrivial bipartitions of n 
parties. There are states that cannot be represented as 
such a mixture. These states, the states that are not 
biseparable, are called genuine miltipartite entangled. 
The problem of recognizing which states are genuine en¬ 
tangled is very important for different applications. It is 
a highly nontrivial task to develop practical conditions 
for detecting genuine multipartite entanglement. In the 
following a single party always corresponds to a single 
electromagnetic mode as known from continuous-variable 
(CV) quantum information theory [Ull]. Many such con¬ 
ditions for n-partite CV systems deal with lower bounds 
for the second-order quantity 

n 

g = (FMf) = = Tr(M7), (1) 

* 

where f = (x, p) is the 2n-vector of position and 



2 


momentum operators and M = is a real, 

symmetric, positive definite 2 n x 2 n matrix and 7 = 
is the covariance matrix of the state, 
({A,i?}) = {1/2){Ab + BA). In fact, a special case 
of matrix M with zero off-diagonal blocks is usually 
used in practice. Denoting X = and P = 

which are symmetric, positive semidefinite 
n X n matrices, the quantity Q in this case reads as 


g = {yAX± + pTpp) = Tr(X7,,) + TiiP-fpp), (2) 


where 72 ,,^ = and jpp = are 

the diagonal blocks of the covariance matrix of the state. 
For example, the classical result of Ref. [3] uses rank-one 

matrices X = ^ and P = = 

gg"’’, where h"'" = (a,a“^) and g"^ = (a, —a“^). The 
works nig use general rank-one matrices X = hh"’’ and 
P = gg"*", where h and g are real n-vectors. In Ref. [ 6 ] 
the quantity Q with 3x3 matrices was used. Refs. [Hi 
deal with the general quantity g. Among the works that 
also use second order moments we mention [nn]- As 
based on second-order moments only, all these criteria 
are sufficient entanglement witnesses for all CV states, 
but necessary and sufficient only for Gaussian states. 

In a realistic setting errors of measurements are un¬ 
avoidable, so to be completely rigorous we need to in¬ 
corporate the possibility of errors into our scheme. We 
do this and formulate a hierarchy of entanglement con¬ 
ditions as convex optimization problems for g in terms 
of the covariance matrix 7 and the information about er¬ 
rors of the measurements. For discrete variables a convex 
optimization approach has been developed in Ref. m- 

The paper is organized as follows. In section Illwe ob¬ 
tain the minimal value of g defined by Eq. p) over all 
quantum states. In Section [lll| we show how to improve 
the lower bound obtained in the preceding section for 
separable states and construct a hierarchy of separability 
conditions in the form of convex optimization problems. 
In Section [Tv] we apply our construction to rank one ma¬ 
trices and demonstrate that it coincides with some pre¬ 
viously known results, so that those results are just a 
special case of our more general approach. In Section [V| 
we give an example of an entangled state with positive 
partial transposition that can be detected by our condi¬ 


tions. Section VI is devoted to an analytical condition 
for genuine multipartite entangled that can be obtained 
from our hierarchy of conditions, in Section [VII| we fur¬ 
ther extend our approach by taking measurement errors 
into consideration and demonstrating that our method 
works in realistic settings as well. Then we give a con¬ 
clusion and provide in appendices all technical details 
missing in the main part of our work. 


II. QUANTUMNESS BOUND 

First of all, we find the minimum of g over all quan¬ 
tum states and then we show how this bound can be 


improved for separable states. To find the minimal value 
of g for a given matrix M we use the Williamson’s theo¬ 
rem [T3]. This theorem states that there is a symplectic 
matrix S such that S'^MS = ( q where A is a diag¬ 
onal n X n matrix. Since each symplectic transform is 
implementable as a unitary transformation |14] . starting 
with a state g we have 


pMi) = i) (3) 

for the appropriately transformed state g'. The minimum 
of g is achieved for g' being the vacuum state, and in this 
case the equality 1/ = TrA takes place. We thus obtain 
that the inequality g = Pv{M^) > TrA is valid for all 
quantum states and it is tight. 

To compute the minimal value of g we need to know 
the diagonal elements Xj of A. These numbers are re¬ 
ferred to as symplectic spectrum of M and they can be 
directly obtained from the matrix M according to the fact 
that ±*Aj are the eigenvalues of JM, where J = ( f) , 
and E is the identity matrix. In our special case of block- 
diagonal M = p) we can get these numbers directly 
in terms of X and P. In fact, we have the equality 

The characteristic equation of this matrix is x(A) = 
det (^x \e) ~ the diagonal blocks commute 

with the off-diagonal ones, according to m this equa¬ 
tion can be simplified as x(A) = det(A^A + XP) = 0. 
Substituting the eigenvalues HXj into this equation we 
see that the diagonal elements Xj satisfy the equation 

det(A:P - X^E) = det{VXPVx - X^E) = 0, (5) 

from which it follows that they are the eigenvalues of 
the symmetric matrix \/ \fXP\fX. The matrices X and 
P can be swapped in this derivation. We arrive to the 
following result: The minimal value of g is given by the 
expression 


min g = Tr\J^/XP^/X = Tr\Jy/PX^/P. ( 6 ) 

§ 

To put it in another way, we have the tight inequality 

Tr(A 7 ,,) +Tr(P 7 pp) > B{X, P) = Tr ^JVxpVx, (7) 

which is valid for all positive semidefinite matrices X and 
P and all correlation matrices "/xx and 7 pp. This inequal¬ 
ity gives some bound on g. The bound has been obtained 
without any assumptions about separability properties of 
the state in question and thus is valid for all multipartite 
quantum states. Our task now is to improve this bound 
for separable states. The improvement will depend on 
the separability kind of the state — the more separable 
groups of modes the state has, the higher is its separa¬ 
bility kind in the partition lattice and the more can this 
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estimation be improved. In Appendix A we also show 
that from the inequality Q one can get a special case of 
Araki-Lieb-Thirring trace inequalities [MIH]. Such in¬ 
equalities play some role in quantum entropy theory, see 

Ref. [H]. 


III. HIERARCHY OF SEPARABILITY 
CONDITIONS 

We first show that pure states with real wave functions 
are enough to minimize Q = Q{X,P). In fact, in terms 
of wave functions the quantity Q reads as 

g = J ((xTAx)|V-(x)| 2 + (VV'*)^P(VV')) dx. (8) 

If we take a general wave function of the form '0(x) = 
where /(x) = |' 0 (x)| is a real wave function 
and (^(x) is the phase, we will see that the first term in 
Eq. ^ does not depend on y)(x), while the other term is 
equal to 


I ((V/)Tp(V/) + dx 

> j{yf)^p{yf)d^. 


( 9 ) 


It follows that for any wave function '0(x) there is a real 
wave function /(x) (the absolute value of the former) 
such that {X, P) >gf (X, P) and thus we can consider 
only pure states with real wave functions. 

If a pure state with real wave function is {/i,..., 7^}- 
separable then 7 pp,ij = (piPj) = 0 if the indices i and j 
belong to different blocks Ii. In fact, separability for pure 
states means that the wave function is factorizable, i.e., 
ip{x) = /(x')g(x"), where, without loss of generality, we 
can assume that x' = (cci,... ,Xk), x" = {xk+i, ■ ■ ■ ,Xn)- 
For 1 < i < k and fc -f 1 < j < n we have 



The fact that '!/'(x) is real is important for the validity 
of the last step. The seeming asymmetry in position and 
momentum operators is superficial — if we worked in mo¬ 
mentum representation and dealt with real wave function 
in that representation we would have (xiXj) = 0 . 

We have just shown that (piPj) = 0 if modes i and j 
are separable, but for the position correlations we can say 
only that they factorize: {xiXj) = {xi){xj). We can get a 
similar conclusion for position moments if we take mini¬ 
mization property into account. If a state \tjj) minimizes 
g then we can also assume that {xiXj} = 0. In fact, tak¬ 
ing the wave function ' 0 o(x) = '(/'(x-l-xo), where Xq = (x) 
is the vector of averages computed for the wave function 
'!/'(x), we get a new wave function with (x)o = 0 and 


go = g — (x)'^X(x) < g. it X is positive definite then 
we must have (x) = 0 and thus {xiXj) = 0 for any separa¬ 
ble state '0(x) that minimizes g. If X is degenerate then 
we can just find a separable state with {xiXj) = 0 that 
minimizes g, but there can be minimizing states that do 
not have this property. Moreover, this new wave func¬ 
tion ■i/'o(x) also has the property that (piPj) = 0 if i and j 
are separable. We see that among {/i,..., /fe}-separable 
states minimizing g we can always find a pure state with 
real wave function for which = lpp[Ii\Ij\ = 0 

for I < i, j < k, i ^ j. The notation A[/|/'], where / and 
/' are sets of indices, is used to denote the submatrix of 
A formed by the intersection of rows with indices in I 
and columns with indices in Using such a state, we 
can improve the lower bound for g. 

Due to the relations (xiXj) = (piPj) = 0 for separable 
i and j we have the equality 


Tr{X'^xx) + Tr(P7pp) = Ti{Xujxx) + Tr(Pv7pp). (H) 


where Xu is obtained from X = {xij)'lj^i by replacing 
its elements corresponding to zero elements of ^xx by ar¬ 
bitrary real numbers Ui subject to the condition that Xu 
is symmetric and positive semidefinite and the same pro¬ 
cedure is applied to P to produce P^ ■ In other words, we 
can replace all elements of the submatrices of the form 
X[Ii\Ij] and P[Ii\Ij], i ^ j, by arbitrary real numbers 
in such a way that the resulting matrices are again sym¬ 
metric and positive semidefinite. This construction is 
better illustrated by an example. For n = 4, consider 
2|134- and l|2|34-factorizable states. In the former case 
we have Ii = {2}, I 2 = {1,3,4}, and the the latter case 
we have Ii = {!}, I 2 = {2}, = {3,4}. The matrices 

Xu corresponding to these two cases are 


fxii Xi3 a:i4\ /xii U2 U 3 \ 


r^cTi X22 [u 2 U3] 


a;i3 

'U2 

2:33 

2:34 


to ^ 

U4 

\a;i4 


2^34 

X44/ 



U5 


X22 U4 U5 


( 12 ) 


In the first case we replace elements of the submatrices 
X[Ji|J 2 ] and X[/ 2 |/i] by arbitrary numbers, and in the 
second case we replace submatrices X[Ji|/ 2 ], X[/ 2 |/i], 
X[/i|/ 3 ], X[/ 3 |/i], X[/ 2 |/ 3 ], X[/ 3 |/ 2 ]. Different subma¬ 
trices are marked by different colors (symmetric parts 
are marked by the same color). The more factorizable 
parts the state has the more elements can be replaced by 
arbitrary numbers. For a completely factorizable state 
we can freely choose all the off-diagonal entries. 

In order not to overload the notation, we fix some kind 
of separability, i.e., some decomposition Z = {/i,..., I/.} 
of the indices, and use it in all the considerations below. 
Applying the inequality Q, we find that for a pure I- 
factorizable state with real wave function, and thus for 
all I-separable states, we have 

g(X, P) > i3x(X, P) = max Tr (13) 

u.v ’ 
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where the optimization is over the points u and v such 
that Xu and Pv are positive definite. For example, for 
a bipartition {/i,/ 2 } with |/i| = m, I/ 2 I = n — m the 
vectors u and v have m{n — m) components, so in this 
case the optimization problem has 2m[n — m) variables. 
The full separability optimization problem has n(n — 1) 
components. 

Note that for the vectors Uq and vq with the corre¬ 
sponding elements of the matrices X and P, respectively, 
we have Xuq = X and Pvq = P and thus for any partition 
Bx{X, P) > B{X, P). This inequality is strict for most of 
the matrices X and P. It follows that the inequality (13) 
gives a stronger bound on Q then does the inequality 
Moreover, in the same way we obtain that if P is finer 
than I then Bx' (X, P) > Bx{X, P), since in this case we 
have more variables to optimize over, and thus the condi¬ 
tion for X'-separable states gives a stronger bound then 
the condition for X-separability. The condition for full 
separability gives the strongest bound of all. We have 
just obtained a hierarchy of separability conditions that 
mirrors the partition lattice in Fig.[^ 

To each kind of separability X corresponds its own 
maximization problem of its own dimension. There are 
2 "“^ — 1 different bipartitions of the indices of an n- 
partite state and many more partitions into three or more 
parts. If for a given state there is a pair of matrices X and 
P such that an inequality (13) is violated, then the state 
is not X-separable. If there are X and P such that the 
inequalities (131 are violated for all bipartitions simulta¬ 
neously, then the state is genuine multipartite entangled. 


IV. RANK-ONE MATRICES 


We now show that the results of R ef. H I are just a 
special case of the general inequalities and give an 
analytical solution of these optimization problems in this 
special case. Consider the rank-one matrices X = 
and P = gg"*". The square root of a rank-one ma¬ 
trix A = aa"^ is given by \fA = aa'^/||a||, so we have 
B{X,P) = |(h, g)|. As a concrete example let us con¬ 
sider four-partite case and X = l|2|34-separable states. 
We are free to change some elements of the matrices X 
and P. Let us just change the sign of the X and P’s 
elements that are marked in Eq. (12|: 


X' = 


/ h\ ±hxh 2 ±/11/13 ±/ii/i4\ 

± hlh2 /I 2 i ^ 2^3 ±/l2^4 
±/ii/i 3 ±/12/13 h^hx 

\ zLhihx i/i2li4 113/14 


hi 


(14) 


For appropriate combinations of signs we can get that 
X' = h'h'"”" and P' = g'g^"^, where the new vectors read 
as h' = (±/ii,±/i 2 ,/i 3 ,/i 4 ) and g' = (±gi, ± 32 , ffa, 54 ), 
and thus 


Pi(X,P) > max|(h',g')| 

= max I ± higi ± h 2 g 2 + hsgs + / 1454 I 
= \higi\ + |/i252| + \h 3 g 3 + hxgxl > B{X,P). 


( 15 ) 


This result can be extended to all n and arbitrary kind 
of separability and coincides with the results obtained in 

Ref. d]. 

Note that the inequality h'^ 7 a; 2 ,h -|- g^'^ppg > |(h,g)| 
is equivalent to the inequalities 


/ Jxx 2 ^®^ 

\2^ ^ppj 


> 0 . 


(16) 


The same idea can be applied to the separability con¬ 
ditions and we arrive to the conditions for X-separable 
states 


f ^xx 2^^\ 
\ 2 'ypp j 


>0, 


(17) 


where Ex is a diagonal matrix with diagonal elements 
equal to ±1 so that the elements with indices in the same 
group Ij have the same sign. These conditions work in 
many cases, but sometimes the more general conditions 
are needed. In the next section we give an example of a 
PPT state that satisfies the inequalities (17) but violates 
the general inequalities (13). 


V. AN EXAMPLE OF PPT STATE 

Consider a four-partite state with the covariance ma¬ 
trix given by [HI [2^ 


'^xx — 


Ipp — 


1 

2 


1 

2 


/2 0 1 0 \ 

0 2 0 -1 

1 0 2 0 ’ 

yo -1 0 2 / 

/I 0 0 -1\ 

0 1-10 
0-140 
V-1 0 0 4 / 


(18) 


The matrix 7 = ^ j is a covariance matrix of a 

quantum state, since 7 -I- (z/2) J > 0. Partial transposi¬ 
tions of (l,3)-kinds (that is transpositions 1|234, 2|134, 
3|124 and 4|123) are negative, and these negativities are 
easily detected by the conditions For example, the 
matrix Pi 1234 reads as 


'E'1|234 — 


/-I 0 0 0\ 
0 10 0 
0 0 10 
V 0 001 / 


(19) 


and the matrix has negative eigenval¬ 

ues. On the other hand, the partial transpositions of 
( 2 , 2 )-kinds are all positive and the matrices ( [T7| ) with 
X = 12|34,13|24,14|23 are positive semidefinite, so the 
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test 0 does not detect any entanglement of a ( 2 , 2 )- 
kind. Nevertheless, the state (18) is entangled in any of 
this kind. 

To demonstrate this, consider the positive semidefinite 
matrices 


X = 


P = 


/ 

X 

0 - 

^/xy 

0 \ 


0 

X 

0 


-^/xy 

0 

y 

0 

V 

0 

y/py 

0 

y / 

1 

V 

0 0 




0 

P y/m 

0 



( 20 ) 


0 

\Vm 0 


0 
9 


where x, y, p and q are some positive numbers to be de¬ 
termined. The quantity G(X, P) — Pt:{X^xx)-\-" hT{P"fpp) 
reads as Q{X, P) = 2x + 2y— 2y/Py -|- p -b 4g — 2^/^. To 
compute Bi 2 \za{X^P)^ note that 

Sl2|34(^, P) > \l (21) 

where the matrices X' and P' are defined to be (since we 
can play with the elements marked by color) 



(x 

O'' 


0 

X , 

X' - 




0 

0 



A 


P' = 


0 y) 


fp 0^0 0 '' 
0 p ,0 0. 
q 0 


flTd 

0 0 


( 22 ) 


0 q) 


and thus for the boundary /Bi 2|34 we have the inequality 
^12|34 > 2(v^+V^)- (23) 

It is easy to check that for the following values of the 
parameters: 


a: = 0.144375, ?/ = 0.084087, 

p = 0.232000, q = 0.039543 


(24) 


we have G{X, P) = 0.435170, while according to Eq. (23) 


we have ^612134 > 0.481359. We see that ^612134 > 
G{X,P), and the PPT state Eq. (18) is 12|34-entangled. 
The construction for the partitions 13|24 and 14|23 is 
similar. 


VI. GENUINE ENTANGLEMENT GONDITION 


The inequalities (13) allow one to test states for en¬ 


tanglement of some kind, but the number of these kinds 
grows extremely fast with the number of parts. We now 
derive an analytical condition for genuine multipartite 
entanglement. It is a single condition that does not re¬ 
quire testing exponentially many bipartitions, however. 


it does not provide the best possible bound, 
the quantity 


Consider 


Gn = 


E 


{{Xi + Xjf+ {p^-Pjf). (25) 


It is the general quantity G with the following matrices 
X = Xn and P = Pn'. 


X„ = 


Pr,. = 


/n 

- 1 

1 

1 

1 \ 


I 

n — 1 

1 

1 


1 

1 

... n — I 

1 

V 

1 

1 

... I 

n — 1/ 

/n 

-1 

-1 

... -I 

-1 \ 


-1 

n — 1 

... -I 

-1 


-1 

-1 

... n—1 

-1 

V 

-1 

-1 

... -I 

n — 1/ 


(26) 


Since matrices and P„ are completely symmetric with 
respect to different parts, it is enough to consider only 
bipartitions of the form 1 ,..., k\k -|- 1 ,..., n for k < n/2. 
We do not change the elements of and in the ma¬ 
trix Pn we set all Vi to 1 (so we change the matrix 
elements from —I to I). Denote the resulting matrix 
by Pnk- The matrices Xn and P^ ^ commute, so that 
TT{\/XnPn is easy to compute 


Tr \/XnPn^k = (n- l)v/n(n- 2) -b 4/3 ^^^ 


where the number is given by 
/?= ^ 


(27) 


(28) 


yj2n — 2 -b y/n — 2 

The minimum of this expression over 1 < fc < n/2 is 
attained for k = 1. We have just obtained the following 
result: Any biseparable state must satisfy the inequality 

4(n — 1) 


Gn > (n-l)\/n(n- 2 )-b- 


(29) 


y^{\/2n — 2 + \/n — 2) 

If this inequality is violated, then the state is genuine 
multipartite entangled. Table summarizes the lower 
bounds of Gn for some n obtained with the analytical 


condition (29) and computed numerically from Eq. (13). 


A similar genuine entanglement condition has been ob¬ 
tained in Ref. [1] in terms of rank-one matrices. The gap 
between quantum bound and biseparability bound there 
decreases as 0 (l/n), where n is the number of r)arts| 21 ). 
In our condition this gap is 0(1) and so it does not tend 
to zero for large n. 


VII. GONVEX OPTIMIZATION AND 
EXPERIMENTAL ERRORS 

To demonstrate entanglement of a kind I we need to 
find a pair of matrices X and P that violate the condition 
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n 

2 

3 

4 

5 

6 

7 

8 

q 

0 

3.46 

8.48 

15.49 

24.49 

35.49 

48.49 

a 

2 

5.00 

10.03 

17.06 

26.07 

37.08 

50.09 

b 

2 

5.46 

10.89 

18.26 

27.59 

38.89 

52.17 

f 

2 

6 

12 

20 

30 

42 

56 


TABLE I. Lower bounds for Qn for different kinds of n-partite 
states. The q row is quantumness bound, (n — l)\/n(n — 2). 
The a row is the biseparability bound given by the analytical 
expression (29l. The b row is the true biseparabili ty bound 
given by the solution of the optimization problem (13l. The 
last row, /, is the full separability bound, n(n — 1). 


(13). In practice, however, we should take into account 


the errors in the measurement of the matrix elements of 
')xx and 7 pp. Assuming that the errors in the individual 
elements of ^xx and 7 pp are independent, the standard 
deviation a{X, P) of G{X, P) is given by the expression 


\x,p)=Y,i 


2 2 
T • <T 

^ XX,IJ 


I 2 2 \ 


(30) 


where <Jxx,ij and (Jpp ij are the standard deviations of 
individual elements of jxx and jpp respectively. So, to 
be on the safe side the right inequality reads as 

£i{X, P) = g(X, P) + s (t(A, P) - Bx{X, P) > 0, (31) 

where s is the level of certainty with which we can claim 
that the state is entangled. We prove that the function 
£x{X,P) is convex on the set of all pairs of semidefinite 
matrices {X,P). 

To prove the convexity of £x{X, P) we have to prove 
the concavity of Bx{X, P) because the convexity of the 
other two terms of £x{X^ P) is obvious. The key element 
of this proof is the fact that Tr \/y'XuPvV^u is jointly 
concave with respect to all four variables X, P, u and v. 
Due to the equality 

Tr \/ VxpVX = min (Tr(A 7 ,,) + Tr(P 7 pp)) (32) 

i'Tpp 


standard argument given in Ref. [22) can be applied here 
to conclude that Bx{X,P) is concave as the maximum 
of a jointly concave function over a convex set. This fin¬ 
ishes the proof of the convexity of the function £x{X, P) 
defined by Eq. (31). 


To violate the inequality (31) we thus have to opti¬ 


mize £x{X, P) over X and P and check whether the op¬ 
timal value is negative or not. This optimization problem 
has n(n -|- 1) variables, the elements of X and P. Since 
the function £x{X, P) is homogeneous, £x{XX, XP) — 
\£x{X,P) for A > 0, it makes sense to put some condi¬ 
tion on the matrices X and P. The simplest is a linear 
condition, for example, the condition Tt{X^xx + P^pp) = 
C, where C is an arbitrary fixed positive constant. We 
thus arrive to the following Z-separability condition: 


min (s cr(A, P) - Bx(X, P) > -C. (34) 


If, for given ^xx, Jpp, ctxx and Upp, this inequality is vi¬ 
olated, i.e., if this minimum drops below —C then the 
state in question is not separable of the corresponding 
kind. If these inequalities can be violated for all bipar¬ 
titions simultaneously (by the same pair of matrices X 
and P), then the state is genuine multipartite entangled. 
The methods to solve convex optimization problems like 
the one given by Eq. (34) are discussed, e.g, in Ref. f^ . 


What s should we choose? Usually, the ’’three-sigma 
rule”, s = 3 is applied [23] • To better understand what 
values of s in Eq. (31) are sufficient to guarantee that our 


results are correct we need to know the probability that 
the result of a measurement lies outside s sigma interval. 
Eor a Gaussian probability distribution with the mean /i 
and the standard deviation a this probability is given by 
the expression 


P(s) = 1 — / ! _ e 2 ,,^ da; = 1 — erf ( —= I , 

K-sa \V2j 

(35) 

which depends only on s. This probability decreases very 
quickly as s growth. The order of values of P(s) for some 
s are shown in the table below. In many cases, the value 


it immediately follows Tr \/y/XP^/X is concave with re¬ 
spect to X and P as the minimum of a family of concave 
(in fact, linear) functions. From the relation 

(0Xi + (1 - 0)X2)'(dUi + (1 - e)u2) 

= 0W((ui) + (l-d)A'(u 2 ), ^ ^ 

0 < 0 < 1, and similar relation for P we derive the joint 
concavity of Tr \/\^XuPvVXu- 

We now have three sets — the set £lx,i x Dp i of points 
(u, v) where Ar((u) and P((v) are positive definite, the 
similar set £lx,2 x Dp _2 for X2 and P2, and the set £lx x Dp 
for A = eXi + (1 - e)X2 and P = ePi + {l- e)P2. In 
general, these are distinct sets, but one can easily see 
that 6£lx,i X Dp^i -I- (I — d)flx ,2 x Dp _2 C D^, x Dp. The 


s 

P{s) 


123456 7 10 

0.32 0.05 10“® 10-* 10“® 10"^ 10"^2 10“^® 


of s = 3 is sufficient (the three-sigma rule). For s > 5 
this probability is negligible, and for s > 6 it is prac¬ 
tically zero. Even if the real probability distribution is 
not perfectly Gaussian it is unlikely to have long tail, so 
Eq. (35) gives a reasonable estimate. Even if this esti¬ 
mate is wrong by several orders of magnitude, provided 
that we have verified violation with s > 6 we are still on 
the safe side. The larger s we set the larger the probabil¬ 
ity of the correct result is, but the more difficult it will 
be to find a violation with such s. From this table we 
can conclude that we should search for a violation with 
s not smaller than 3 and not larger than 6 — the event 
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of getting the right result outside of six sigma interval is 
practically impossible. 

We see that if we can violate our inequalities with s = 6 
then it practically guarantees that our conclusion is cor¬ 
rect. In Appendix B we demonstrate the usefulness of 
our approach by applying to some four-, six- and ten- 
partite realistic states. We also demonstrate that the 
four-partite state in question is genuine multipartite en¬ 
tangled. 


VIII. CONCLUSION 

We have developed a method to test continuous- 
variable multipartite states for arbitrary kinds of entan¬ 
glement. Our approach allows both numerical and an¬ 
alytical treatment. Numerically, it reduces to a convex 
optimization problem, which allows fast and accurate so¬ 
lution. We have shown that it is very efficient at detect¬ 
ing ordinary entanglement and can detect genuine mul¬ 
tipartite entanglement in a reasonable amount of time. 
Analytically, it allows to reproduce (and thus generalize) 
some known results as well as to obtain an analytical 
genuine multipartite entanglement condition. With our 
approach we can easily obtain a trace-class inequality, 
which is difficult to get in a direct way. 


Appendix A: Trace inequalities 


If matrices X and P commute then the right-hand side 
of Eq. ([^ reduces to Ti'/X'/P = Tr '/P^/X. This ex¬ 
pression, Tt{VxVp), is a lower bound for Q independent 
of commutation properties of X and P. For a real wave 
function /(x) we have the equality 

g{X,P)= y(u'^(x)Xu(x)-I-v'’’(x)Pv(x))dx, (Al) 


where the vector fields u(x) and v(x) are defined via 
u(x) = /(x)x and v(x) = V/(x). We can write this 
equality in a more compact form as 


g{X,P) = y(||u(x)f+ ||v(x)f)dx, (A2) 

where the new vector fields are defined via u = \/Xu 
and V = y/Pv. Now we can estimate g as follows: 

g{X,P)>2\ y(u(x),v(x))dx| 


= 2 

From the relation 


^VA/P(V/)(x)/(x)dx 


(A3) 


J Uj{x.)vk{x)dx = J dx = (A4) 


we get the inequality g{X,P) > Tt{'/X'/P). 


We thus have two lower bounds for g — the tight one 
is given by the inequality Q and the other one, not nec¬ 
essarily tight, have just been obtained with the help of 
Cauchy-Schwarz inequality. Since the tight bound is the 
best bound possible, we derive the following inequality 
for a pair of positive definite matrices X and P: 

Tr \I'/XP'/X > Tr(VA/P) = Tr( v^/P</A). 

,(A5) 

This inequality is a special case of Araki-Lieb-Thirring 
trace inequalities |16H18j . which also have quantum me¬ 
chanical background and read as 

Tr(AiAsAiA)’'9 > Tr(A’'A5M’'A)9, (A6) 


where A and B are arbitrary positive definite matrices, 
q > 0 and 0 < r < 1. The c ase of <7=1 and r = 1/2 
corresponds to the inequality (A5). 


Appendix B: Application to realistic matrices 


Ixx - 


It happens that rank-one matrices work surprisingly 
well. Consider the four-partite state that was analyzed 
in Ref. [5]. It has the following covariance matrix: 


/ 1.09921 0.16092 -0.17609 -0.84831\ 

0.16092 0.40938 -0.16060 -0.18963 

-0.17609 -0.16060 0.46060 0.04319 

V-0.84831 -0.18963 0.04319 1.06419 / 

/1.09921 0.35533 0.36439 0.91386\ 
0.35533 0.92282 0.57440 0.43388 
0.36439 0.57440 1.04339 0.34868 
VO.91386 0.43388 0.34868 1.06419/ 


Ipp — 


(Bl) 


^ XT, - 


The standard deviation matrix reads as 

/0.00327 0.01041 0.00894 0.00647\ 

0.01041 0.00822 0.01848 0.01899 
0.00894 0.01848 0.00861 0.01345 
VO.00647 0.01899 0.01345 0.00549/ 

/0.00458 0.01009 0.02767 0.04289\ 

0.01009 0.01023 0.02101 0.02085 
0.02767 0.02101 0.01466 0.01955 
V0.04289 0.02085 0.01955 0.00455/ 

A very simple way to search for violation of the condi¬ 


(B2) 


tion (31) is to randomly generate 4-vectors h and g and 
check whether this condition is violated by the rank-one 
matrices X = hh"”" and P = gg"*", and if it is, how strong 
the violation is. Then just record the maximal observed 
violation. As a measure of violation we use the quantity 


s = 


Bi{x,p)-g{x,p) 

a{X,P) 


(B3) 

















This approach requires only simple matrix algebra ma¬ 
nipulations, which can be done very efficiently with tools 
like Intel Math Kernel Library. A simple parallel Fortran 
program has been written and run on a low-end 4-core 
desktop PC. The total time to test all seven possible bi¬ 
partitions in this four-partite case is 4 minutes (using all 
four cores available). Table [Tl| compares our results with 
those obtained in Ref. [S] . We see that for the state under 
study our approach is superior to that of Ref. [5] (which 
uses a genetic algorithm to find the best violation), since 
it is simpler and gives better results, though, as we have 
mentioned before, from practical point of view all viola¬ 
tions larger than 6 are of the same value. 


Bipart. 

Violation 

h 

Ref. [8] 

4m 

g 

1|234 

20.93 

26.48 

(1.97,-0.01,0.49,1.88) 

(1.14,0.18,-0.20,-1.03) 

2|134 

13.17 

18.08 

(-0.40,-1.99,-1.48,-0.74) 

(-0.15,-1.80,1.24,0.68) 

3|124 

11.21 

16.10 

(0.31,1.93,1.62,0.32) 
(0.23,1.65,-1.46, -0.19) 

4|123 

21.06 

26.57 

(1.77,0.18,0.44,1.82) 

(-0.96,-0.18,-0.12,1.21) 

12|34 

24.34 

27.79 

(1.91,0.42,0.75,1.90) 

(-1.06,-0.60,0.50,1.19) 

13|24 

23.52 

26.17 

(-1.52,-0.17,-0.54,-1.58) 
(0.78,-0.09,0.24, -0.97) 

14|23 

4.66 

9.72 

(0.22,1.65,0.18,0.77) 
(0.28,1.34, -0.37,-0.95) 


TABLE II. The comparison of the violation of the separability 
condition for the state given by the covariance matrix (Bl| 


with measurements errors given by Eq. (B21. The second col¬ 
umn shows the results obtained in Ref. [8J, the third column 
lists the results obtained by randomly testing the inequalities 
(311. The total time to perform all seven tests is 4 minutes. 


We now apply our technique to the six-partite state 
also considered in Ref. [5]. We have performed two runs 
of our program on the same hardware as in the previous 
case, one with a smaller number of random trials and 
the other with 200 times more trials. The first run takes 
approximately 4 minutes to perform all 31 tests, the other 
one takes around 12 hours. As Table. [m] demonstrates, 
in this case the optimization based on a genetic algorithm 
gives somewhat larger violations. On the other hand, we 
do not know what computational resources were used to 
perform that optimization and how much time it took. 
As we have already said, all violation larger than 6 are of 
the same practical value and our method produced much 
better violations in just a few minutes on a low-end PC. 

The last state considered in Ref. [S] is a ten-partite 
state. It has been reported that the smallest violation 
of 1.1 was obtained for the bipartition 1,10|23456789. 
The corresponding probability to get wrong result is 
P(l.l) « 0.27, and it is not small enough to conclude that 
the state under study is not 1,10|23456789-separable. 
Randomly generating vectors h and g, we have found 
that the inequality (31) for this kind of separability can 


Bipartition 

Violation 

Ref. ig 

12h 

4m 

1123456 

40.086 

40.111 

37.500 

2| 13456 

36.185 

34.097 

33.967 

3| 12456 

20.274 

19.715 

17.683 

4| 12356 

20.010 

18.869 

16.953 

5| 12346 

27.146 

26.680 

22.979 

6| 12345 

49.220 

48.077 

44.187 

12|3456 

53.541 

53.085 

50.142 

13|2456 

45.569 

45.958 

41.684 

14|2356 

44.789 

42.163 

42.509 

15|2346 

45.282 

40.410 

38.565 

16|2345 

31.177 

29.995 

25.686 

23|1456 

40.158 

38.636 

37.462 

24|1356 

37.698 

36.633 

31.716 

25|1346 

35.256 

31.106 

29.877 

2611345 

47.016 

42.269 

40.199 

34|1256 

24.833 

22.592 

20.694 

35|1246 

28.794 

25.927 

23.390 

36|1245 

50.193 

48.021 

44.561 

45|1236 

30.629 

28.950 

27.092 

46|1235 

51.500 

50.153 

47.510 

56|1234 

56.080 

55.390 

51.521 

123|456 

56.661 

55.474 

54.345 

124|356 

54.402 

52.666 

52.044 

125|346 

50.653 

48.953 

47.993 

126|345 

28.957 

26.470 

22.905 

134|256 

47.675 

46.956 

46.154 

135|246 

47.279 

43.109 

40.016 

136|245 

34.237 

29.400 

26.527 

145|236 

47.331 

42.957 

41.199 

146|235 

35.340 

34.072 

30.970 

156|234 

39.487 

38.426 

36.880 


TABLE III. The comparison of the violation of the separabil¬ 
ity condition for the six-partite state considered in Ref. [8]. 
Provided that the separability condition can be so strongly vi¬ 
olated for all bipartitions it is absolutely unnecessary to test 
other kinds of separability, i.e. kinds with partitions of modes 
into three or more gronps. 


be violated with s = 3.65. The corresponding probability 
P(3.65) < dTO”"* is much smaller and gives a strong con¬ 
fidence that the state is 1,10|23456789-entangled. The 
vectors are 

h = ( - 0.65, -1.22, -0.21,0.01,0.365, 

- 0.32,0.25,0.28,-0.88,-0.94), , ^ 

g = (0.22,-1.04,-0.12,-0.04,0.42, ^ ^ 

- 0.34,0.22,0.30,-0.56,0.73). 

The violations of other kinds of biseparability are all 
larger than 3, so the standard three-sigma test is passed 
for all bipartitions. The violations reported in Ref. [S] 
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show some strange behavior — the violation for full sep¬ 
arability is smaller than violations of some more coarse 
kinds. But this may be an artifact of an implementation 
of the genetic optimization algorithm. 

Up to now it has been shown that the four-partite state 
with the covariance matrix (Bl) is not separable for any 
fixed kind of separability. We demonstrate that this state 
is genuine entangled. To do this we need to find a pair 
of matrices X and P that simultaneously violate the in¬ 
equalities (31) for all bipartitions. First, we have tried to 
violate these inequalities with rank-one matrices. It took 
nearly one day, but we were able to find a pair of vectors 


the maximum is attained at 


/ 0.39234 

( - 0.10873 0.158136 0.116524\ 

[-0.1087J 
0.158136 
\[ 0.116524 

0.88526 0.09450 0.09080 

0.09450 0.58391 0.20795 

0.09080 0.20795 0.39504 / 

/ 0.22992 (- 0.11914 0.113758 0.083761]\ 

-0.119H 

0.113758 
'1 0.083761 

0.52598 -0.32316 -0.16699 

-0.32316 0.39949 0.06971 

-0.16699 0.06971 0.31242 / 


h = (0.31, -1.93, -0.17, -0.18) 
g = (0.30,1.48,-0.57,-0.53), 


(B5) 


that violate the conditions (31) for all seven biparti¬ 
tions, and the minimal violation is 3.15 (for the bipar¬ 
tition 1|234). The corresponding probability P(3.15) = 
1.6 • 10“^ is relatively small to conclude that the state 
under study is genuine entangled. 

The approach with a simpler conditions works but it 
takes a lot of time and it just marginally passes the three- 
sigma test. Using the general matrices we can do better. 
The sketch of our approach is as follows. We use a vari¬ 
ant of the steepest gradient method. According to this 
method, to optimize a convex function one has to go in 
the direction opposite to the gradient of the function. 
Here we have several functions to be optimized at once, 
and each has its own gradient. We start by generating a 
pair of random positive definite matrices X and P and 
compute the gradients of all seven target functions. If 
the directions of these gradients are not strongly scat¬ 
tered then we can take the average of the gradients, go 
in the opposite direction and still improving all our func¬ 
tions simultaneously. If the gradients point into nearly 
opposite directions then we cannot proceed this way, so 
we stop and generate a new random pair of matrices. We 
do this until we find proper matrices X and P or give up 
after some prescribed number of attempts. Following this 
approach, in a few hours we found the following pair of 
matrices for the four-partite state with the covariance 
matrix (Bl |: 


/ 0.39234 

-0.20267 

0.24691 0.30527\ 


-0.20267 

0.88526 

0.09450 0.09080 


0.24691 

0.09450 

0.58391 0.20795 


\ 0.30527 

0.09080 

0.20795 0.39504yi 


/ 0.22992 

-0.13140 

-0.00477 -0.11723\ 

-0.13140 

0.52598 

-0.32316 -0.16699 

-0.00477 

-0.32316 

0.39949 0.06971 

V-0.11723 

-0.16699 

0.06971 0.31242 / 


For these matrices we have G = 1.47484 and a = 0.01947. 
The bound Bx{X,P) for different bipartitions is pre¬ 
sented below. The elements of the matrices that were 
optimized over are highlighted. For the bipartition 1|234 


and is equal to Hi| 234 (A, P) = 1.65474. For the biparti¬ 
tion 21134 at 


A' = 


0.24691 
V 0.30527 


-0.01993 


P' = 


-0.00477 
V-0.11723 


- 0.0143? 
0.02340 


) 0.24691 

0.30527 \ 

[ - 0.03586 

-0.01993] 

0.58391 

0.20795 

0.20795 

0.39504 / 

) -0.00477 

-0.11723\ 

[ - 0.01432 

0.02340] 

0.39949 

0.06971 

0.06971 

0.31242 / 


and is equal to H 2 |i 34 (A, P) = 1.66193. For the biparti¬ 
tion 31124 at 




P' = 


/ 0.39234 -0.20267 [0.22149 ]0.30527\ 
-0.20267 0.88526 -0.01671 0.09080 


(0.22149 -0.01671) 0.58391 [0.24154] 
V 0.30527 0.09080 [0.24154] 0.39504/ 

/ 0.22992 -0.13140 fO.02836 
-0.13140 0.52598 


■0.07629 


-0.11723\ 
-0.16699 

[0.02836 -0.07629] 0.39949 [0.12842] 
\-0.11723 -0.16699 [0.12842] 0.31242 / 


and is equal to P 3 |i 24 (A, P) = 1.56935. For the biparti¬ 
tion 4|123 at 


A' = 


P' = 


0.15483\ 

0.04047 

0.23966 


/ 0.39234 -0.20267 0.24691 

-0.20267 0.88526 0.09450 

0.24691 0.09450 0.58391 

\ [0.15483 0.04047 0.23%(r) 0.39504/ 

/ 0.22992 -0.13140 -0.00477 [0.05094 1\ 
-0.13140 0.52598 -0.32316 -0.05608 

-0.00477 -0.32316 0.39949 [o. 12953 J 
V [ 0.05094 -0.05608 0.1295T) 0.31242 / 


and is equal to P 4 |i 23 (A, P) = 1.63974. For the biparti- 
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tion 12134 at 


tion 14123 at 


X' = 


P' = 


/ 0.39234 -0.20267 
-0.20267 0.88526 


0.19766 -0.022^ 

V [0.11649 0. 03156 


a 


/ 0.22992 -0.13140 
-0.13140 0.52598 


f0.11949 -0.02001] 
0.06522 0.02362 


[0.19766 

0.11649\ 

-0.02260 

0.03156] 

0.58391 

0.20795 

0.20795 

0.39504/ 

0.11949 

0.06522^ 

-0.02001 

0.02362] 

0.39949 

0.06971 

0.06971 

0.31242/ 


and is equal to Bi 2 \ 3 i{X,P) = 1.81056. For the biparti¬ 
tion 13124 at 




P' = 


/ 0.39234 [ - 0.10013] 0.24691 [0.12997\ 

( -0.10013] 0.88526 [ - 0.03695) 0.09080 
0.24691 f - 0.036951 0.58391 f 0.25436] 

V [0.12997] 0.09080 [0.25436] 0.39504/ 

/ 0.22992 [ - 0.07273] -0.00477 [ 0.06225 ] \ 
[ -0.07273] 0.52598 [ - 0.05209] -0.16699 
-0.00477 [ - 0.05209] 0.39949 [0.17734] 

\ [0.06225] -0.16699 [0.17734] 0.31242 


and is equal to Bi 3 \ 2 i{X,P) = 1.74993. For the biparti- 


/ 0.39234 

[ - 0.11435 

0.18571] 


0.30527 \ 

r- 0.11431 

0.88526 

0.09450 


- 0.02361) 

0.18571, 

0.09450 

0.58391 


0 . 25307 ] 

V 0.30527 

[ - 0.02360 

0.25307] 


0.39504 / 

/ 0.22992 

[ - 0.09531 

0.05497] 

-0.11723 

[ - 0.0953: 

. 0.52598 

-0.32316 

- 0.0217] 

0.05497 
^_ 

-0.32316 

0.39949 


0 . 12643 ] 

\ -0.11723 [ - 0.02171 

0.12643] 

0.31242 


and is equal to Bii\ 23 {X,P) = 1.56114. The smallest 
number among these maximums is the last one, 1.56114, 
so we have 


Biix,p)-g{x,p) 

a{X,P) 


>so = 4.43199 


for all bipartitions I simultaneously. The corresponding 
probability is P(so) < 10“®, which is almost two orders 
of magnitude smaller than for the vectors h and g we 
found before, so one can be pretty sure that the state 
under study is genuine entangled. 
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